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ABSTRACT 

We present a versatile family of model galactic outflows including non-uniform mass and energy 
source distributions, a gravitational potential from an extended mass source, and radiative losses. 
The model easily produces steady-state wind solutions for a range of mass-loading factors, energy¬ 
loading factors, galaxy mass and galaxy radius. We find that, with radiative losses included, highly 
mass-loaded winds must be driven at high central temperatures, whereas low mass-loaded winds can 
be driven at low temperatures just above the peak of the cooling curve, meaning radiative losses can 
drastically affect the wind solution even for low mass-loading factors. By including radiative losses, 
we are able to show that subsonic flows can be ignored as a possible mechanism for expelling mass and 
energy from a galaxy compared to the more efficient transonic solutions. Specihcally, the transonic 
solutions with low mass-loading and high energy-loading are the most efficient. Our model also 
produces low-temperature, high-velocity winds that could explain the prevalence of low-temperature 
material in observed outflows. Finally, we show that our model, unlike the well-known Chevalier & 
Clegg (1985) model, can reproduce the observed linear relationship between wind X-ray luminosity 
and star formation rate (SFR) over a large range of SFR from 1 — 1000 M 0 /yr assuming the wind 
mass-loading factor is higher for low-mass, and hence, low-SFR galaxies. We also constrain the allowed 
mass-loading factors that can fit the observed X-ray luminosity vs. SFR trend, further suggesting an 
inverse relationship between mass-loading and SFR as explored in advanced numerical simulations. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — galaxies: star formation 
— radiation: dynamics — X-rays: galaxies 


1 . INTRODUCTION 

Due to the large variety of physical phenomena occur¬ 
ring on a vast range of scales, galaxy formation is very 
challenging to model theoretically, and the full complex¬ 
ity of many processes still needs to be unraveled. Cur¬ 
rently, many models of galaxy formation without some 
form of feedback form stars too early and too quickly, 
consequently overestimating the number of low-mass 
galaxies and their stellar ages compared to observations 
(White & Frenk (1991), Oppenheimer & Dave (2006), 
Piontek & Steinmetz (2011), Hopkins et al. (2012a), Vo- 
gelsberger et al. (2013)). One possible solution is to in¬ 
troduce a more efficient galactic wind feedback model in 
these systems. Galactic winds from supernova rich star- 
burst regions and active galactic nuclei transfer mass and 
energy away from the regions of wind development and 
enrich the intergalactic medium with metals, eventually 
suppressing star formation in the galaxies. 

There is a wealth of observational evidence for galactic 
winds (Pettini et al. (2003), Shapley et al. (2003), Martin 
(2005), Veilleux et al. (2005), Steidel et al. (2010), Gen- 
zel et al. (2011)). However, large-scale galactic outflows, 
both their intrinsic properties and their effects on galaxy 
formation and evolution, are poorly understood. Sim¬ 
ulations of thermally driven winds (Chevalier & Clegg 
(1985), Ciotti et al. (1991), Li et al. (2015)), radiation 
driven winds (Murray et al. (2005), Murray et al. (2011)), 
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and cosmic ray driven winds (Everett et al. (2008), Ev¬ 
erett et al. (2010), Uhlig et al. (2012), Hanasz et al. 
(2013), Salem & Bryan (2014)) have been developed and 
successfully reproduce observations of galaxies, includ¬ 
ing the Milky Way and the well-known starburst galaxy, 
M82. However, these physical processes are at too small 
scales to be properly accounted for in cosmological sim¬ 
ulations. Instead, hydrodynamical simulations of galaxy 
formation (Schaye et al. (2010), Stinson et al. (2012), 
Vogelsberger et al. (2012), Aumer et al. (2013), Hop¬ 
kins et al. (2013), Christensen et al. (2015), Muratov 
et al. (2015), Schaye et al. (2015)) typically capture these 
feedback processes using sub-resolution physics guided 
by parameterizations from observations. These outflow 
prescriptions cannot be specific to one or a few galax¬ 
ies; instead, they must be general enough to encompass 
galaxies on many scales. It is logical, then, to propose 
simplified but versatile models for galactic winds and at¬ 
tempt, with those models, to determine what the most 
important parameters are for driving a wind and how 
these simplified models compare to observations. 

In Section 2, we describe our complete re-working 
of the Chevalier & Clegg (1985) model, including non- 
uniform mass and energy source distributions, a gravita¬ 
tional potential from an extended mass distribution, and 
radiative losses. In Section 3, we solve for a set of crit¬ 
ical points at which the flow transitions from subsonic 
to supersonic. We do this after scaling the system of 
steady-state wind equations by GM/R, as shown in Sec¬ 
tion 3.1, and using the logical mass-loading and energy¬ 
loading parameterizations of Section 3.2. We show the 
critical points for a number of different mass-loading and 
energy-loading factors in Section 3.3. In Section 3.4, we 
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describe the effects of radiative losses on the wind solu¬ 
tions. We find that radiative losses can be important for 
all ranges of mass-loading factors, even the least mass- 
loaded winds, which can generally be driven at lower 
temperatures near the peak of the cooling curve. We 
also show in Section 4 that radiative losses can severely 
decrease the efficiency in which winds expel mass from 
a galaxy. For subsonic winds, as shown in Section 4.2, 
radiative losses sharply decrease the temperature, and 
hence energy per mass expelled in the wind, at a radius 
dependent on the energy-loading factor. For transonic 
solutions, we show in Section 4.1 that the least mass- 
loaded outflows with the highest energy-loading factors 
are the most efficient. In Section 4.3, we also note the ex¬ 
istence of very low temperature, high velocity transonic 
outflows in our model due to radiative losses. Finally, 
we compare, in Section 5 the resulting relationship be¬ 
tween outflow X-ray luminosity and star formation rate 
(SFR) to recent observations. Unlike the classical Cheva¬ 
lier and Clegg model, which has been shown to produce 
a quadratic relationship between luminosity and star for¬ 
mation rate (Zhang et al. (2014)), including gravity and 
supposing higher mass-loading fractions for lower star 
formation rates can fairly well produce the observed lin¬ 
ear relationship. Specifically, we show that, to match 
the linear behavior of the observations, the allowed mass¬ 
loading factors generally must decrease monotonically as 
star formation rate, and hence galaxy mass, increases. 


2 . MODEL 

One of the standard models for galactic winds was pro¬ 
posed by Chevalier & Clegg (1985), hereafter referred to 
as the CC model. The model supposes a spherically sym¬ 
metric wind, does not include a gravitational potential, 

and supposes a constant mass injection q = ^ and con¬ 
stant energy injection Q = ^ out to a radius R. For 
r > R, q = Q = 0. For a spherical geometry, the mass 
conservation, momentum, and energy equations are 
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where u is the wind velocity. For the CC model, q{r) = 
qo = constant for r < R, and d^/dr = 0. This model 
was originally proposed to describe M82, a well-studied 
starburst galaxy with an extremely strong wind. X-ray 
fits indicate high terminal wind speeds of up to ~ 1400 — 
2200 km/s (Strickland & Heckman (2009)), so far above 
the escape velocity of the galaxy that there is no need to 
include gravity in a model to describe M82; hence, the 
CC model omits it. 

Parker (1958) showed that there are infinitely many 
wind solutions, all corresponding to different boundary 
conditions. Many of these are “breeze” solutions where 
the wind accelerates up until the “Parker point” and then 
decelerates to zero velocity again, always staying sub¬ 
sonic. We discuss subsonic outflows in Section 4.2, but 


we show that radiative losses inhibit their efficiency, lim¬ 
iting the mass flux to certain radii before a sharp drop¬ 
off in temperature. The most effective type of wind for 
transferring matter out of a galaxy is the transonic so¬ 
lution that monotonically increases in velocity, becomes 
supersonic at a critical radius r^, and has an asymptotic 
velocity considerably greater than the subsonic solutions. 
For now, we just consider transonic solutions. 

One can solve the CC mass conservation, momentum, 
and energy equations for the Mach number of this tran¬ 
sonic wind as a function of radius. The radius, R, at 
which there is a discontinuity in mass and energy input 
is also the critical radius, Tc for the CC model, i.e. it is 
the radius at which the wind switches continuously from 
subsonic at r < i? to supersonic at r > R. A notice¬ 
able, non-physical consequence of this model is that the 
velocity derivative is infinite at radius R. This is a conse¬ 
quence of the jump in mass and energy injection at that 
radius. 

To instead get an analytic velocity solution at r = R, 
we modify the mass injection to now be a function of r, 
namely 

q{r) = go(l - ^) (4) 

out to a radius R. We choose this idealized mass injec¬ 
tion profile instead of, e.g. an exponential or power law 
profile, because it allows for tractable, analytic calcula¬ 
tions while still retaining the logical and observed behav¬ 
ior that more mass per volume should be injected to the 
wind near the center of the galaxy. The result of using 
a non-uniform mass source instead of the CC constant 
mass source is shown in Fig. 1 for a polytropic wind 
with polytropic index 7 = 5/3 in the absence of gravity 
and energy injection. Note that including a non-uniform 
mass source shifts the critical point so it is no longer 
exactly at r = R. This will prove to be important, as 
there can actually be infinitely many critical points cor¬ 
responding to solutions with different central conditions. 
For each polytropic wind, such as that shown in Fig. 1, 
the shifted critical point for our q(r) profile is always very 
close to r/R = 1 regardless of the factor K in the wind 
equation of state P = Kp^. In addition, in the region 
r > R where there is no mass or energy addition, the 
wind expands exactly as it does for the CC model. One 
can see in Fig. 1 that the asymptotic behavior for each 
model is very similar. Therefore, we do not believe that 
the results in this paper are sensitively dependent on the 
q(r) profile we have chosen. 

It should also be noted that the CC model energy equa¬ 
tion, eqn. (3), is not the same as the one we use in our 
own model, though it is very similar. In the CC model, 
each term on the left hand side has units of energy per 
time per volume, hence it is logical to evolve those quan¬ 
tities with an energy per time per volume factor, Qcc- 

We choose instead to determine the temperature evo¬ 
lution by the first law of thermodynamics; dU = dQ + 
Pdp/p^, where dQ — TdS, and Q is not to be confused 
with the CC source term in their energy equation. Using 
P = (pf — l)pU and assuming a steady state with radial 
flow u, the temperature equation is 


ks dT 
m dr 
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Figure 1. The CC model is not analytic at the critical radius R. 
This problem stems from having a constant mass source qo. Using 

2 

a mass source such as q(r) = qo{l — gives a finite slope at the 
now shifted critical point. 


and the pressure is described by 


dP 2dp . {l-l)pdQ 

where = 'ykBT/m is the sound speed squared, and 
we set 7 = 5/3 in all cases studied in this paper. If 
dQ/dt = 0, i.e. we have no additional energy injection, 
then the resulting flow is adiabatic. The CC model, on 
the other hand, requires this extra energy injection to 
drive a wind because the central pressure would be zero 
without energy addition, whereas our model always has 
a non-zero pressure at the galactic center. 

A logical way to shape the additional energy addition 

term, dQ/dt, is to give it a 1 — profile, as is the case 
for the mass injection. Therefore, we define an energy 
per mass parameter rj such that 


P 


dQ 

dt 


= VQ 


(7) 


Equation (1) is solved in the Appendix, giving 


P = 


r{q) 

3u ’ 


( 8 ) 


where (q) is the mean q inside radius r. 

We also extend the model to include an extended mass 
distribution chosen for simplicity to have constant mass 
density 


Pm = 3 M/{A-kB?^) for r <Rm (9) 

While this is not fully realistic, having the star-forming 
mass more centrally concentrated than the total mass 
replicates some features of central starbursts. Inclusion 
of a non-uniform mass density with an extended halo is 
the subject of current work. 

Our gravitational potential is then defined in two re¬ 


gions: 


GM , 

2Rm^Rl ^ 

for r < Rm 

( 10 ) 

GM 

$ =- 

for r > Rm 

( 11 ) 


In this paper, we always set Rm = R, i-e. we tie the 
end of the mass distribution to the end of the mass and 
energy source distribution; however, the sonic point is 
not restricted to this radius, as we will see in Section 3.3. 

3. METHODOLOGY 
3.1. Scaling by 

It is useful to scale the main wind equations by the 
quantity which represents the magnitude of the 

gravitational potential energy per unit mass at r = i?. 
By doing so, each of our solutions to the scaled equations 
will in fact give us a multitude of solutions for various 
gravitational potentials. In addition, we write our equa¬ 
tions in terms of the dimensionless radius f = r/R (see 
eqns. (A 8 ) and (A9) for the momentum and sound speed 
equations, respectively, in terms of scaled variables). The 
resulting scaling for the main wind quantities is then 


/ i?3 

( 12 ) 

^ \ GM 

c,=CsVgWR 

(13) 

u = u\/ GM/R 

(14) 


where the tilde denotes the scaled variable. The energy 
per mass addition also scales accordingly: ry = fjGM/R. 

3.2. Parameterizations 

To include mass-loading in our modified CC model, we 
parameterize the mass-loading efficiency by 

M = ^SFR(Mo/yr) (15) 

Then our mass-loading per volume factor, qo, is calcu¬ 
lated such that M = f^qo(l — r^/R‘^)dV. 

qo = = 1.60 X 10-3^/3SFR (16) 

^ttR3 

fixing R = 200 pc for each galaxy, regardless of the 
galaxy’s mass and SFR. Then our gravitational potential 
is calculated as GM/R. We consider galaxies of dynam¬ 
ical mass 10 ® — lO^^M©, and the assumption that this 
entire mass is confined to a radius of 200 pc from the 
center is quite unreasonable for many galaxies. We are, 
however, considering only starburst galaxies, for which 
10 ^ — 1 O®M 0 of molecular gas can be present near the 
galactic center. In future simulations, we will change the 
galaxy radius and see what effect that has on our wind 
solutions. We would also like to study the effect of a 
more realistic mass distribution, e.g. including a dark 
matter halo, etc. Including even a simple gravitational 
potential, which is no worse than including a point mass 
potential, allows us to model the many galaxies that do 
not have a very strong, M82-like wind. For these galaxies. 
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the inward pull of gravity is quite important, regardless 
of whether the gravitational potential is strictly realistic. 

We parameterize energy addition as 

E = a' Esn (17) 

Esn = ei^SFR is the energy per time injected into 
the wind from supernovae where e = lO^^esi ergs is 
the energy injected by a single supernovae, and v = 
*^ioo/(lOOM 0 ) is the number of supernovae per unit 
mass of star formation. We assume ^loo = 1, mean¬ 
ing that one supernovae occurs for every IOOMq of stars 
produced. For a Salpeter initial mass function (IMF), 
= 1-18, and for a Chabrier IMF, ^loo = 1.74 (Lei- 
therer et al. (1999), Strickland & Heckman (2009)). Then 
a = a'esi^^ioo is how we parametrize the energy injection, 
following Zhang et al. (2014) to most easily compare to 
their results using the CC model. In the CC model, this 
means Q = a'EsN/V, whereas in our model, this pa¬ 
rameterization means the energy per mass injected to 
the wind is 

r] = 5.032 X lO^^ai/ioo (18) 

and 77 = r]/{GM/R). 

The fj = 0 case corresponds to a polytropic wind, i.e. 
a wind obeying a law of the form T oc P (X . For 

this scenario, supernovae heat the gas to a certain cen¬ 
tral temperature, Tq; however, after this central heating 
occurs, the wind is then driven solely by adiabatic ex¬ 
pansion and advection, not from any additional energy 
deposition into the wind following a 1 — r^jP? profile, as 
would be the case for 7) > 0 in our model. 


3.3. Determining the Critical Points 

To get the wind properties (velocity, temperature, den¬ 
sity, etc.) as a function of radius, one must simultane¬ 
ously solve the mass continuity, momentum, and energy 
equations, i.e. eqns. (1), (2), and (5). Due to the non¬ 
linearity of the equations, solving for the steady-state 
solution is best done numerically by integrating inwards 
to r « 0 from a known critical point, where the critical 
velocity, temperature, etc. are known, and also outwards 
to greater radii. 

Finding the critical point for a polytropic wind is 
straightforward. When adding additional energy (77 > 
0 ), it is not as straightforward to find the critical 
point. We follow the method outlined in Tamers & 
Cassinelli (1999): One writes the momentum equation 
as F{r,u,T, du/dr,dT/dr) = 0 and solves the following 
set of equations 


F{r,u,T, du/dr, dT/df) = 0 (19) 


dF dF dh OF df 
dr du df dT df 


( 20 ) 


where du/dr can be found using L’Hospital’s rule at the 
critical point. Because F = 0 all along the solution curve 
by definition, 

dF 

dF/df = G+ —u" = 0 (21) 

dw 


all along the solution curve, as well. One can solve for 


the second derivative of velocity: 


At the critical point. 


G{f,u,f) 

OF 

du' 


dF 

dxf 



(22) 


(23) 


So for u” to be well-defined at the critical point, 
G{f,u,T) should be exactly zero at the critical point. 
If this is true, then du/df is continuous at the critical 
point, meaning that u{r) will be a smooth function at 
the critical point. We would like this to occur, so we 
search for critical points by solving F = 0 and G = 0 
simultaneously. 

Solving eqns. (19) and (20) without specifying the ini¬ 
tial conditions for the desired transonic solution yields 
the same overlapping solution curves, i.e. an infinite 
number of smooth critical points, each corresponding to 
different initial conditions, as seen in Fig. 2. Fig. 3 gives 
the set of critical points for varying values of 77 , showing 
that the critical velocity generally shifts upwards with 
an increase in energy addition. It should be noted that 
the critical point locations for non-radiative winds do not 

depend on our mass injection parameter qo; however, the 
2 

^ ~ W is important. 



r/R 


Figure 2. M = IO^Mq, R = 200 pc. Many polytropic (7 = 5/3), 
transonic wind solutions, each corresponding to a different cen¬ 
tral sound speed and, hence, different central temperature. Cen¬ 
tral wind temperatures range from 1.1 x 10^ — 1.3 x 10^ K. Each 
wind is driven purely thermally, with no additional energy or mass 
injection. Wind velocities monotonically trend higher for higher 
central wind temperatures; therefore, solutions with high critical 
and asymptotic velocities can only be achieved with a high central 
temperature. 


3 . 4 . Effect of Radiative Losses 

We further examine two polytropic, non-radiative wind 
solutions (shown in Fig. 4), and we show the ratio of the 
cooling timescale, Gooi = ikbT /riA, to the dynamical 
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F = 0, G = 0 Solution Curve 



Figure 3. Critical point solution sets for varying energy per mass 
inputs rj ranging from 0 to 10. The critical points are found by 
solving the equations F = 0 and G = 0. These two curves are 
the same; therefore, there are many critical points with smooth 
velocity profiles, each for a different central wind temperature. To 
get unsealed velocity, multiply scaled velocity by yjGM/K. The 
right axis shows the corresponding central wind temperature for a 
galaxy of lO^M© and rj = 0.0 


timescale, tdyn = Jq dr/u, of the solutions for various 
star formation rates and mass-loading efficiencies. This 
provides a first order estimate as to whether radiative 
losses should be included. The cooling timescale trends 
lower for lower SFR and also for higher /3; however, we 
also see that radiative cooling is not nearly as important 
if we use a scaled solution with higher initial tempera¬ 
ture (solution 2) than if we use one with lower initial 
temperature (solution 1). Nevertheless, it is important 
to include radiative losses for all temperature ranges. 

With radiative losses included, our mass continuity 
equation stays the same, but our momentum and tem¬ 
perature equations now contain radiative loss terms: 


pu 


du 

dr 


= -qu - c. 


,dp 

dr 


5d> ,.F^A(T) 

■ + (7- 1)- 


dr 


m^u 


(24) 


ks dT 
m dr 
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dQ 

dt 


+ (7- 1) 


ksTu dp 
Trip dr 


( 7 - 1 ) 


pA(T) 

17772 


(25) 


where p/Tfi = n is the number density, to = 1.021 x 10“^^ 
g is the mean mass of a wind particle, and A(r) is the 
optically thin radiative loss function. 

Instead of using a more detailed cooling curve such 
as from Schure et al. (2009), we choose for simplicity 
to use an analytic fit very similar to that of Imada & 
Zweibel (2012), given by eqn. (26), which fits the peak 
near 5 x 10® K fairly well and also captures the behavior 
for both higher and lower temperatures, assuming colli- 
sional ionization equilibrium (CIE) 


A(r) = Aol0®(®) (26) 


where d(T) = logiQ(T/ro)j To = 2x10® K, and we choose 


Aq = 1.5 X 10 ergs s ^emr’ 

0(6») = 0.40 - 3 

-f 6.2/[exp(1.50 -F 0.08) -F exp(-(0 -F 0.08))] (27) 

A plot of the analytic solution vs more detailed calcula¬ 
tions from Schure et al. (2009) is given in Fig. 5. 

For now, we just consider transonic solutions for ra¬ 
diative winds. The system of steady-state, radiative 
wind equations is given by eqns. (24) and (25). As 
for the non-radiative case, we scale these equations by 
GM/R and solve for the critical points, i.e. we solve 
F{f,u,f,du/df,df/dP) = 0andG = + + 

= 0. As for the non-radiative case, we find that 
F = 0 and G = 0 define the same solution curves, which 
are given in Fig. 6 for various mass-loading factors jd. 
Looking at Fig. 6, we also see that the critical velocities 
tend greater with increasing /3; an increased critical ve¬ 
locity is only achieved with an increased initial temper¬ 
ature, assuming the initial velocity always approaches 
u = 0 at r = 0. Therefore, we see that more heav¬ 
ily mass-loaded winds (higher jd) require a higher initial 
temperature to achieve any transonic solution. This is 
more explicitly shown in Fig. 7. 

Recall from Fig. 4 that greater mass-loading for the 
same initial temperature (same solution curve) implies 
a more radiative wind due to an increase in density; 
however, because the cooling time is proportional to the 
temperature, higher temperature winds are generally less 
radiative. Therefore, since high-/3, radiative winds can 
only become transonic if given a high temperature, the 
radiative losses incurred from heavily mass-loading the 
winds are somewhat lessened by the higher temperature. 
Conversely, since it is possible to have low-temperature, 
low-/? transonic solutions, radiative effects can also be 
important in low mass-loaded winds. 

We show this in Fig. 8 by comparing two transonic 
winds with the same initial temperature: one with ra¬ 
diative loss physics included and one without for various 
mass-loading factors. The resulting velocity is decreased 
by about a factor of two at large radii, the tempera¬ 
ture decreases by a factor of 2 — 4, and the density is 
slightly increased for large radii after starting off with a 
lower central density for the radiative solutions. A plot 
of temperature is given in Fig. 8. 

4. EFFICIENCY OF WINDS IN EXPELLING MASS 

We have modified the widely-used Chevalier and Clegg 
wind model by including an extended mass source term, 
a simplified gravitational potential, and radiative losses 
using an analytic cooling curve that fairly well encapsu¬ 
lates the behavior of more detailed cooling curves at all 
temperatures. With these logical modifications, we be¬ 
lieve our model includes all the necessary physics while 
still being quite simple, allowing us to probe the fun¬ 
damental questions of importance to observers and the 
galaxy formation community: What conditions are nec¬ 
essary for a thermally driven wind? How efficiently can 
mass be expelled from a galaxy by a thermally driven 
galactic wind? Are subsonic winds or transonic winds 
more efficient outflows? 

The energy per unit time driven outward by a wind is 
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Two Solutions Analyzed in Further Detail Cooling Timescale vs Dynamical Timescale 




Figure 4. Two transonic solutions {left panel) and the ratio of the cooling timescale to the dynamical timescale of each solution over a 
range of SFRs for various 0 and M = IO^Mq {right panel). The lower velocity and lower sound speed solution is hereafter called solution 
1, and the higher velocity and sound speed solution is hereafter called solution 2. To get unsealed velocity, multiply scaled velocity by 
yjGMjK. The cooling timescale for winds in our model trends lower for higher SFR and also for higher /3, both of which increase the 
density of the wind; however, we also see that radiative cooling is not nearly as important if we use a scaled solution with higher initial 
temperature (solution 2) than if we use one with lower initial temperature (solution 1). 



Figure 5. To self-consistently include radiative losses in our 
model, we choose to use the analytic fit from Imada & Zweibel 
(2012), which fairly well fits the peak near 5 x 10^ K as well as the 
structure at higher temperatures. It also captures the steep drop 
in emissivity for temperatures below 10^ K; however, our assump¬ 
tion of CIE should not be trusted below about 5000 K. For low 
temperatures, one would need to include a non-equilibrium model. 
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Figure 6. SFR = IMq/ yr, galaxy mass M = IO^Mq. These 
are solution curves for critical points of radiative winds. One can 
observe that higher critical velocities, and hence higher initial tem¬ 
peratures, are required to achieve a transonic solution with higher 
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given by 


f 

= M 
Air 



7 P 


7 - 1 p 


(28) 


excluding gravitational potential and radiative loss terms 
that decline with distance and do not contribute asymp¬ 
totically. 

Using M = /3SFR, 


47r 


^SFR 


2 ^ 



We define the energy per mass quantity 


e = 


L 9 

-.u + 



(29) 


(30) 


This effectively quantifies how efficiently the winds can 
expel mass and energy from the galaxy. 


4.1. Efficiency of Transonic Outflows 

We have shown that for each mass-loading factor fl 
and energy factor a, there are infinitely many transonic 
solutions, each for a different central wind temperature, 
which can drive mass to large radii. The question we 
would like to ask is “for the same central wind tempera¬ 
ture, what are the trends in a and fl that lead to the most 
efficient mass expulsion from a galaxy?” We look at this 
by considering a set central temperature of 5 x 10^ K, a 
SFR of lOOM 0 /yr, and a galaxy mass of IO^Mq, which 
would be fairly typical for an active starburst galaxy, 
and we vary a and fl. We then calculate e at a radius 
of r/R — 3, or 600 pc, at which point the kinetic en¬ 
ergy dominates the energy flux, for each radiative wind 
solution to track the energy per mass of the wind. The 
solutions considered and a contour plot of energy per 
mass for the /3 < 0.4 are shown in Fig. 9. 

At high /3, radiative losses are more important due to 
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Figure 7. SFiJ = IMQ/yr, galaxy mass M = 10®Mq. Wind 
temperature profiles corresponding to a representative sample of 
critical points in Fig. 6 with no additional energy input (o; = 0). 
In general, the low-/3 transonic solutions can be driven at lower 
central temperatures, whereas the high-/? solutions can only be 
driven at higher temperatures. The existence of low temperature 
winds for which the wind temperature plummets to close to 10^ K 
will be discussed in Section 4.3. 


Radiation’s Effect on Temperature 



Figure 8. Plot of temperature for /? = 0.04,0.4,4. As shown 
in Fig. 6, a more mass-loaded (higher /?) wind requires a higher 
central temperature to achieve a transonic solution. Comparing 
radiative to non-radiative winds, we see that the effect of radiation 
is more significant in higher mass-loaded winds (due to higher den¬ 
sity in the wind) than the lower mass-loaded winds; however, the 
difference is not extreme because the high-/? winds are also those 
with high temperature. Conversely, the low-/? solution has lower 
density; however, it also has a lower temperature, meaning that 
radiative effects can still be important even if the wind is not very 
dense. 

the wind being denser, resulting in a lower temperature 
and lower velocity wind. Therefore, e decreases as /3 
increases. This can be offset by a large a, however, as 
e generally increases as a increases. We also see that, 
unless the wind has a high a > 0.5, solutions with (3 > 
0.4 drop in temperature below 10^ K very sharply within 
600 pc, at which point our assumption of CIE becomes 
invalid, and our calculation of e shouldn’t be trusted. 
Including radiative heating at temperatures 10^ K will 


be left to further work. 

Something to note about the sound speed profiles 
in Fig. 9 is that, for high a, the sound speed, and 
hence temperature, increases near the galactic center. 
This is because energy is being injected at a high rate, 
and advection and radiative losses cannot cool the gas 
fast enough. This is an interesting result because the 
l — profile for energy injection would seem logical, 

as you expect more supernovae to occur near the galactic 
center. We hope to include heat conduction in a future 
model, which should decrease this rise in temperature. 

4.2. Ejficiency of Subsonic Outflows 

We have so far only considered transonic solutions, but 
there are a wealth of subsonic solutions that could also 
expel mass from a galaxy, albeit at a slower pace since the 
asymptotic velocity is much lower. In Fig. 10, we show 
a set of subsonic solutions with various a and the same 
central temperature, (3 — 0.04, SFR = lM 0 /yr, M = 
IO^^Mq. The lowest sound speed curve corresponds to 
the lowest velocity curve, for which a = 0. For greater a, 
the temperature and velocity curves are higher. Because 
radiative losses are greatest at low temperatures near 
the cooling curve peak, the low a (lower temperature) 
curves are most affected by radiative losses, causing a 
sharp drop in temperature at some radius. For the a = 0 

solution, this sharp drop occurs at roughly r/R = 1.8, 
2 

meaning that + ~ ^ drops by orders of magnitude 

at « 360 pc. In Fig. 10, we roughly track the point at 
which temperature and, hence, e drops sharply for winds 
of various a. We see that winds with higher a can expel 
mass to greater radii, but only so far. The transonic 
solution, however, is the only solution that can efficiently 
expel mass to very large radii. Therefore, we conclude 
that transonic solutions are indeed the most important, 
especially those with high a and low /3. 

4.3. Low Temperature, High Velocity Outflows 

We notice from Fig. 6 that the sonic points for an out¬ 
flow can be closer to the galactic center, but still near the 
edge, depending on the wind’s central temperature. This 
is because low temperature winds cool more rapidly than 
high temperature winds further away from the cooling 
curve peak. Fig. 11 shows a transonic solution that cools 
very rapidly to temperatures in the range of 10 ^ — 10 ^ 
K, at which point the wind goes over the peak of the 
cooling curve, and the temperature profile flattens con¬ 
siderably. We have encountered many other such winds 
in our models that drop within a few hundred parsecs to 
even lower temperatures, at which point the assumption 
of CIE breaks down. These low temperature winds seem 
especially common in starburst galaxies with high SFR, 
such as in Fig. 9, where even winds of /3 « 0.5 can drop 
to temperatures < 10^ K very quickly unless the wind 
is also highly energy-loaded with a > 0.5. Despite our 
assumption of CIE breaking down at such low tempera¬ 
tures, we believe that radiative cooling in our model can 
explain the prevalence of cool, fast outflows from galax¬ 
ies. This is in agreement with Silich et al. (2011), as well 
as the recent work of Thompson et al. (2015) who see the 
same drastic radiative cooling for high-/3 winds using a 
different modification of the CC model. For low (3, our 
results also qualitatively agree in that the wind expands 
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Figure 9. Plot of various radiative wind solutions (left panel) and their efficiencies [right panel) according to our measure, e, of energy 
per mass driven outward by the wind. Each wind has a central temperature of 5 x 10^ K, a SFR of lOOM 0 /yr, a galaxy mass of IO^Mq, 
and a galaxy radius oi R = 200 pc. (3 ranges from 0.04 to 4, and a. ranges from 0 to 0.5. Winds with the highest velocity (blue curve) 
also have the highest sound speed (red curves). These high-velocity, high-temperature winds are the most efficient, and correspond to the 
upper left corner of the efficiency plot in the right panel. At high /3, radiative losses are more important due to the wind being denser, 
resulting in a lower temperature and lower velocity wind. Therefore, e decreases as j3 increases. This can be offset by a large a, however, as 
e generally increases as a increases. For winds with /3 > 0.4 and a < 0.5, the temperature drops below 10“^ K within 600 pc of the galactic 
center, making our assumption of CIE invalid and any calculation of e meaningless. A few of these solutions are plotted with dashed lines 
in the left panel. 


To = 1.5 X 10^ (K), M = 




Figure 10. M = lO^^M©. Set of a few subsonic solutions with the same central temperature T = 1.5 x 10^ K, /3 = 0.04, but with different 
a [left panel) ranging from 0 to 0.05. Also a plot of wind extent, i.e. how far these winds penetrate outside a galaxy until the temperature 
(and hence e) drops considerably [right panel). The sound speed curves (green) and velocity curves (blue) in the left panel monotonically 
increase with increasing a. 


approximately adiabatically, as radiative cooling is not 
as important (see Thompson et al. (2015), Fig. 2). They 
also find that photoionization equilibrium (PIE) may be 
quite important in the temperature range of 10^ - 10^ K, 
at which an upturn in A(T) occurs, and photoionization 
heating can balance cooling. The gas can then expand 
at approximately constant temperature and continue to 
radiate (also see Silich et al. (2003), Silich et al. (2004)). 
This, of course, is only possible if f3 is low enough to 
even produce a transonic solution instead of radiating too 
strongly and producing a galactic fountain. The maxi¬ 


mum j3 that results in a transonic wind depends on the 
SFR, and a range of (3 that can produce plausible winds 
matching the observed X-ray luminosity vs SFR relation¬ 
ship will be discussed in Sec. 5. Of course, winds at tem¬ 
peratures below at which point photoionization 

becomes important, do not radiate in X-rays. Thus, our 
treatment of gas at this low temperature has no impact 
on the agreement we find between the observed and cal¬ 
culated relationship between X-ray luminosity and SFR. 
It may be important for other emission diagnostics. 
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Figure 11. An example of a low-temperature wind with /3 = 0.04, a = 0 from an M = IO^Mq, SFR = lM 0 /yr galaxy. Velocity and 
sound speed are plotted {left panel) as well as temperature (right panel). The wind temperature decreases very rapidly as the temperature 
nears the peak of the cooling curve, reaching the sonic point at just rdR = 0.7725, or 154.5 pc assuming R = 200 pc. Once the temperature 
falls below the peak, the emissivity drops considerably very quickly, hence the quick flattening of the temperature profile and the quick 
deceieration of the wind. Radiative losses then become negligible, and the wind propagates adiabatically. 


5. COMPARING TO X-RAY LUMINOSITY 
OBSERVATIONS 

Finally, we would like to test our model against recent 
observations and try to deduce a general relationship be¬ 
tween mass-loading factor and SFR (or galaxy mass). 
Following the work of Zhang et ah (2014), we calculate 
the theoretical X-ray luminosities, predicted by the 
CC model and predicted by our more general, modified 
model for a wide range of star formation rates. X-ray 
luminosities of actual star-forming galaxies are observed 
to have a linear correlation between and SFR (Mineo 
et ah (2014)); however, Zhang et al. (2014) showed that 
the CC model gives a quadratic relationship oc SFR^. 
This is because the CC solution predicts that n cx SFR, 
and oc n^. In addition, for a CC-like wind to have 
a plausible X-ray luminosity, the mass-loading factor 
P = should be less than unity (Zhang et al. (2014)). 
Otherwise, the predicted X-ray luminosity is orders of 
magnitude greater than observed values. This is in con¬ 
tradiction to many galaxy formation codes, which require 
a mass-loading factor between /3 = 1 — 10 to reproduce 
the observed galaxy stellar mass function for low-mass 
galaxies. For example. Bower et al. (2012) compare many 
models of wind feedback, typically finding that high /3 
values give the best fit to the mass function; Barai et al. 
(2013) assume a constant /3 = 2 in their simulations; 
and Puchwein & Springel (2013) find that ,5 = 2 for a 
wind velocity of 484 km/s reproduces the low-mass end 
of the galaxy stellar mass function quite well assuming 
that wind velocity decreases and mass-loading increases 
with decreasing galaxy mass. 

The most recent observationally determined relation¬ 
ship between the total X-ray luminosity and SFR in star¬ 
forming galaxies is given by Mineo et al. (2014). In the 
0.5 - 8 keV band. 


« (4.0±0.4)x lO 3 ®ergs-V(M 02 /r-i) (31) 


Let’s first consider a wind without energy input {a = 
0). We consider the same two transonic solutions from 
Fig. 4 and re-scale the scaled velocity, sound speed, 
and density functions. We get the temperature T = 
mc^/(fykb) where kb is Boltzmann’s constant, 7 = 5/3, 
and m = fj-mn is the mean mass per wind particle as¬ 
suming the mean molecular weight /i = 0.61 for solar 
abundances. 

The X-ray emission from our wind solutions is then 

i:^ai- 8 kev) = / (32) 

The electron and hydrogen number densities, rtf, and 
uh, respectively, are obtained from the wind density 
by n{r) = p{r)/m and then nn = xn = 0.71n 
for solar abundances. We use the XSPEC package 
(https://heasarc.gsfc.nasa.gov/xanadu/xspec/ ; Version 
12.8.2) to calculate the cooling function Z) 

of a hot plasma in GIF, and in keeping with the XSPEC 
documentation, we assume = 1.2n^f. This ionization 
level is actually a function of temperature, but it is al¬ 
most exactly 1.2 for all temperatures that contribute to 
emission in the 0.5-8 keV range. 

Only some fraction, /d < 1 of this diffuse emission will 
be from the hot wind. In M82, fg « 0.1 (Strickland 
& Heckman (2009)); however, many of the galaxies in¬ 
cluded in Mineo’s sample are not M82-like. In fact, as 
shown in our results in Fig. 12, including a linear rela¬ 
tionship between galaxy mass and SFR, as sketched in 
Fig. 1 of Mineo et al. (2014), the predicted Lx — SFR 
relationship is no longer strictly quadratic as was pre¬ 
dicted by the CC model. The linear galaxy mass to SFR 
relationship we used in Fig. 12 was 

M = 10^° X 0.1 X SFR(M 0 ) (33) 

which is a rough fit to the sample in Mineo et al. (2014). 
It should be noted that the Mineo et al. (2014) sample 
excludes many starburst galaxies that can have low mass 
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but high SFRs and may lie on a relationship similar to 
eqn. (33) but with a different constant of proportional¬ 
ity. We have found that assuming a linear relationship 
between SFR and galaxy mass is key to reproducing the 
X-ray luminosity vs SFR relationship observed by Mineo 
et al. (2014). 

Even without energy input (a = 0), Fig. 12 shows a 
linear correlation between and SFR for various ranges 
of SFR depending on which transonic solution we choose 
from Fig. 4. Using a fairly low central wind tempera¬ 
ture, for low SFR galaxies, temperatures at many radii 
are too low to be picked up in the 0.5 - 8 keV band. Con¬ 
sequently, the low SFR region exhibits very low X-ray 
luminosities. If we choose a solution with higher central 
temperature, the low SFR galaxies produce a more linear 
relationship. If the solution curve has a central temper¬ 
ature great enough, the high SFR galaxies produce wind 
temperatures outside the 8 keV limit, causing the total 
X-ray luminosity in that band to fall off rapidly. There¬ 
fore, the galactic central conditions can play a large role 
in what the — SFR correlation is in this 0.5 - 8 keV 
band. 

As seen in Fig. 12, if we use a higher mass-loading fac¬ 
tor /3, the X-ray luminosity shifts upwards for each SFR 
because the wind density increases; however, radiative 
losses can also be important for all /?. The red stars 
show the effect of including radiative losses to points 
along the /3 = 0.4, solution 2 and the (3 = 0.04, solu¬ 
tion 1 curves. It should be noted, however, that in the 
range of SFR « 20 — lOOMQ/yr, the central temper¬ 
atures of the non-radiative /3 = 0.04, solution 1 curve 
cannot drive a radiative wind with /3 = 0.04. At this 
range of SFR, M = /3SFR is high enough, even with a 
low /3, that higher temperatures are required to drive a 
radiative wind. Therefore, radiative wind luminosities 
are not compared to non-radiative wind luminosities in 
that SFR range. 

Generally, radiative losses increase the X-ray luminos¬ 
ity of the wind in the temperature range of 10® — 10^ K. 
In this range, the emissivity is greater for lower temper¬ 
atures. As seen in Fig. 8 , including radiative losses gives 
the wind an overall lower temperature than for the non- 
radiative wind, meaning the emissivity can be factors of 
ten greater depending on how significant the radiative 
losses are. Also, as the wind propagates outwards from 
r « 0 , some thermal energy that would have otherwise 
been converted to kinetic energy is now lost to radia¬ 
tion. Consequently, the wind doesn’t spread outward as 
quickly and, instead, clumps up, resulting in an increase 
in density near the galactic center. This increase in den¬ 
sity and decrease in temperature both serve to enhance 
the radiative loss, and hence increase the X-ray luminos¬ 
ity. The resulting increase is fairly marginal, though, and 
doesn’t change the linear — SFR behavior. 

It has been observed (Chisholm et al. (2014)) that 
higher mass galaxies (and hence higher SFR galaxies 
according to the Mineo et al. (2014) trend) have lower 
mass-loading factors. This has also been seen in hydro- 
dynamical simulations of galactic winds (Hopkins et al. 
(2012b)). Using j3 = 0.4 for low-mass (SFR < 2 OM 0 /yr) 
galaxies and (3 = 0.04 for high-mass (SFR > 2 OM 0 /yr) 
galaxies, coupled with using a solution curve with higher 
scaled initial temperature for low-mass galaxies, we can 


produce a linear correlation for almost the entire range 
of SFRs (see Fig. 13). 

It should be stressed that we have shown just one such 
combination of mass-loading factor, SFR, and wind solu¬ 
tion that produces a linear — SFR relationship. Many 
such combinations exist. It is important, though, that 
one cannot produce a linear — SFR relationship across 
a large range of SFR without switching the wind solution 
and/or /3 for various ranges of SFR. 

In fact, we find that only certain values of /3 are al¬ 
lowed if we would like to construct a linear — SFR 
relationship from radiative outflows while keeping the 
central galactic temperature to a reasonable value. In 
Fig. 14, we look at the space of mass-loading values that 
can fit the observed X-ray luminosity trend scaled by 
fd = 0.01,0.1, andl.O. For each outflow, as shown in Fig. 
7, the central temperature must be above a certain value 
dependent on the mass-loading factor for that wind. This 
sets a lower bound on [3 because if a wind is not mass- 
loaded enough, it can only radiate brightly if the temper¬ 
ature is near the peak of the cooling curve. These low 
central temperatures cannot always be achieved, though. 
Similarly, we can set an upper bound on j3. For highly 
mass-loaded winds, the temperature must be high, i.e. 
away from the cooling curve peak, to limit the luminosity. 
At some point, we consider this central galactic temper¬ 
ature requirement to be unreasonably high. We choose 
this temperature to be = 5 x 10® AT, so only winds 
with central temperatures below this value are consid¬ 
ered. 

For higher SFR, more mass is being injected into the 
wind for the same f3] consequently, as seen in Fig. 14, 
there must be a decreasing trend between Pmin and /3max 
versus SFR, and hence versus galaxy mass, to limit the 
luminosity and keep the central temperature under the 
5x 10® K limit. As shown in Puchwein & Springel (2013), 
increasing /? for low-mass galaxies helps reproduce the 
observed galaxy stellar mass function in their simula¬ 
tions. Therefore, we believe that the positive effect this 
scaling has on fitting the expected — SFR relationship 
provides further validation, along with observations, to 
use this assumption in future galaxy formation simula¬ 
tions. 

For higher luminosities (for greater fd), the range of 
j3 shifts upward, as a denser wind will radiate more 
strongly; however, it is interesting that, for a high fd = 
1.0 and a low fd = 0-01, the range of [3 gets squeezed at 
high SFR until, at SFRs > 250, the target luminosities 
can no longer be achieved by our model for any value 
of p. The fd = 0.1 trend, however, can be reproduced 
for the entire range of SFR, albeit with the range of /? 
again being squeezed at high SFRs. This squeezing oc¬ 
curs because, for high SFR, the luminosity trends higher 
and higher; therefore, if P is too low, the wind will not 
radiate to the necessary extent. This means the slope 
of the minimum p curve must become less negative, all 
while the maximum P curve, which is constrained by the 
5 X 10® K limit on central wind temperature, continues 
to decrease with roughly the same slope. 
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Figure 12. Plot showing the wind X-ray luminosity vs SFR for our model and the observed trend of Mineo et al. (2014) scaled by 
(blue lines). Results show that increasing the mass-loading factor /3 shifts the luminosity upwards at all SFRs, and varying the temperature 
of the wind (solution 1 has lower temperature; solution 2 has higher temperature) shifts the region in which the Lx — SFR relationship 
is linear. Using a fairly low central wind temperature, for low SFR galaxies, temperatures at many radii are too low to be picked up in 
the 0.5 - 8 keV band. Consequently, the low SFR region exhibits very low X-ray luminosities. If we choose a solution with higher central 
temperature, the low SFR galaxies produce a more linear relationship. If the solution curve has a central temperature great enough, the 
high SFR galaxies produce wind temperatures greater than 8 keV, causing the total X-ray luminosity in that band to fall off rapidly. The 
red stars show the effect of including radiative losses to points along the (3 = 0.4, solution I and the /3 = 0.04, solution 2 curves. The result 
is a marginal increase in luminosity. 
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Figure 13. Using solution 1 and /3 = 0.04 for SFR > 20, solution 2 and ^ = 0.4 for SFR < 20. By increasing the mass-loading factor 
and initial scaled wind temperature for low-mass galaxies, we can roughly create a linear relationship between X-ray luminosity and SFR. 
Since we do not have an exact value for for each galaxy in the sample, we cannot use this plot to constrain /3; however, the possible 
relationship that scaled velocity and mass-loading efficiency should be higher for low-mass galaxies gains some validation from this plot in 
that the expected linear correlation between X-ray luminosity and SFR is at least reproduced. As shown in Puchwein ^ Springel (2013), 
this increase in /3 for low-mass galaxies can help reproduce the observed galaxy stellar mass function in simulations. 
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Figure 14. Assuming = 0.01, 0.1, and 1.0 and a maximum 
central wind temperature of 5 x 10^ K, these are the minimum 
and maximum 0 that can reproduce the target X-ray luminosity 
at each SFR following the trend in the 0.5-8 keV range from Mineo 
et al. (2014) scaled by fd. Note that for SFR > 250, the fd = 0.01 
and fd = 1.0 luminosities cannot be reproduced by our model. 
Overall, there is a clear inverse relationship between allowed 0 and 
SFR, and hence between 0 and galaxy mass, as we assume a direct 
relationship between SFR and galaxy mass. 


6 . CONCLUSIONS 

We have re-worked the Chevalier & Clegg (1985) model 
to include non-uniform mass and energy source distribu¬ 
tions, the gravitational potential for an extended mass 
distribution, and radiative losses. After scaling our 
steady-state equations, we generate many different tran¬ 
sonic wind solutions, each for a different central wind 
temperature, for an infinite number of galactic masses 
and radii. Using this set-up, we can easily explore the 
space of central temperature, mass-loading factors, and 
energy-loading factors. 

In Fig. 4, we show the ratio of the cooling timescale 
to the dynamical timescale of our wind model for various 
star formation rates and mass-loading efficiencies. This 
provides a first order estimate as to whether radiative 
losses should be included. The cooling timescale trends 
lower for lower SFR and also for higher /3; however, we 
also see that radiative cooling is not as important if we 
use a scaled solution with higher initial temperature. 

When we include radiative loss physics, we find that 
generating a transonic solution is more expensive for 
more heavily mass-loaded winds, i.e. the wind must 
have a higher central temperature. For these high - 
/3 transonic winds, then, the cooling time is generally 
greater, hence the radiative losses incurred by increasing 
density are somewhat lessened by the higher tempera¬ 
ture required to get the transonic solution in the first 
place. Conversely, low-/3 winds do not require as much 
thermal energy to drive them, i.e. they can be driven 
at lower temperatures closer to the peak of the cooling 
curve. Therefore, radiative losses can be important for 
low mass-loaded winds, as well, despite having a lower 
density. When radiative losses are significant, we find 
that including radiation actually increases the X-ray lu¬ 
minosity compared to a non-radiative wind with the same 


central temperature. 

With radiative losses included in our model, we also 
consider the efficiency of various transonic solutions, 
namely which a and /3 parameters give rise to transonic 
solutions with the greatest energy flux per mass out of 
the galaxy. We find that, for high /3, radiative losses 
are very important and decrease the asymptotic temper¬ 
ature and velocity of the wind. However, efficiency can 
be recovered if more energy is injected. 

Also, given the infinite number of possible subsonic 
solutions, we consider how efficient subsonic winds can 
be in expelling mass from a galaxy. We find that, due 
to radiative losses again, these winds lose their energy 
per mass very quickly as the wind cools towards tem¬ 
peratures near the peak of the cooling curve. We de¬ 
fine the “wind extent” as how far the subsonic solution 
propagates before this large drop in e, and we find that 
winds with higher energy injection can extend further 
than other subsonic solutions but still only so far. Tran¬ 
sonic solutions, on the other hand, can extend to infinite 
distances in our model, albeit with a density that de¬ 
creases as 1/r^ outside the galaxy radius. Therefore, we 
conclude that transonic solutions, especially those with 
low (3 and high a, are the most efficient in driving mass 
from a galaxy. These low-/3 radiative solutions are also 
the least expensive to drive in that they have the lowest 
central temperature requirement. 

We can also generate many transonic wind solutions 
that decrease in temperature very rapidly due to radia¬ 
tive cooling, possibly explaining the existence of cool, 
fast outflows from galaxies. One example for a relatively 
low /3 = 0.04 is given in Fig. 11. 

We also compare the predicted X-ray luminosity of var¬ 
ious outflows to recent 0.5 — 8 keV observations for a wide 
range of SFRs. This provides a first observational test 
for our model compared to the CC model. We find that, 
when the wind temperature at all radii gives energies be¬ 
tween 0.5 - 8 keV, i.e. when the central wind conditions 
are right to do so, our modified CC model with a non- 
uniform mass source and extended mass distribution pro¬ 
duces a linear correlation between total X-ray luminosity 
and SFR for various ranges of SFR. Choosing a higher 
central temperature and higher mass-loading factor for 
low-mass galaxies, we can produce this linear correla¬ 
tion over the entire range of SFR from 1 — lO^M 0 /yr. 
We present only one combination of mass-loading and 
wind solution resulting in a linear — SFR relation¬ 
ship; however, many combinations can achieve this, as 
long as mass-loading and solution are not fixed over the 
entire range of SFR. Specifically, we find that, assuming 
various scaling factors fd, there is a range of /3 for each 
SFR that can fit the observed luminosity. The allowed (3 
values typically decrease as SFR, and hence galaxy mass, 
increases. This is supported by recent observations that 
find an inverse relationship between mass-loading and 
galaxy mass (Chisholm et al. (2014)), and it lends sup¬ 
port to many galaxy formation simulations which utilize 
higher mass-loading factors for low-mass galaxies to help 
reproduce the observed stellar mass function at the low- 
mass end. 
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APPENDIX 


Then, u = and our continuity equation. 


which was ^ ^ = q{r)^ becomes 


1 d{pur‘^) 

Qf 


= q{f)y/R^/GM 


(A6) 


_ <q>r _ <q>r 




3 u 


Here, p = 

The pressure equation is similarly changed. ^ = 

2 ^ , ( 7 -i)p ^ becomes 
dr ^ U dt Liecomes 


dr ® dr 


(7-1) R GM 


GM V R 


-m 


(A7) 


Using this scaling, our main non-radiative equation 
with scaled variables becomes 


The solution to the continuity equation with spherical 
symmetry, eqn. (1), is 


which can also be written 
P = 

where 


r 0 


3 

4 ) 

- / s q{s)ds 
^ Jo 

(Al) 

— qu — 

r{q) 

3m ’ 

(A2) 

Similarly, the scaled 

del _ 3(7 - l)7i7g 

/ s‘^q{s)ds 

Jo 

(A3) 

dr < q> 1 

, 3(7-1)c2m 

r-2 

1 ~ ~ 

< q > r 



<q> 

3m 




^ (7-1) 


fjq (A8) 


—2<q> < q > r du 


3m 


3m2 dr 


(A9) 


is the mean q inside radius r. 

The outward force is provided by the pressure gradient 
(here let’s just consider thermal gas pressure). For an 
ideal gas equation of state 

dP _ ksT dp ^ pkB dT 
dr m dr m dr 

The temperature is determined by the first law of ther¬ 
modynamics; dU = dQ + Pdp/p^ (where dQ = TdS and 
Q is not to be confused with the CC source term in their 
energy equation). Using P = (7 — l)pU and assuming a 
steady state with radial flow m, the temperature equation 
is given by eqn. (5), and the pressure described by eqn. 
(6)- 

The main equation that is integrated numerically is the 
steady state radial momentum equation, eqn. (2). 

Substituting the expression for dp/dr and plugging in 
the result of the continuity equation, our main velocity 
equation, without radiative losses included, is 


The process of scaling by is similar when radiative 
losses are included; however, since the cooling depends 
on temperature, we cannot completely scale out of 
the radiative loss terms, meaning that we must specify 
the quantity before we generate a wind solution. For 
non-radiative winds, could be cancelled out of each 
term in the scaled momentum equation, allowing us to 
generate a scaled solution and then specify after the 
fact. 


du cl 

P'^TTG — 2) 
dr 


-qu{l + 



( 7 -I) ^ 

M ^ dt 
9$ 2cl < q > 
dr 3m 


(A5) 


It is logical to scale the above equation by We 

also scale r by R, i.e. f = r/R, q{f) = qo{l — and 

< q{r) >= go(l - 3/5f^). 

To add energy sources to our model, we define an en¬ 
ergy per mass function ij such that P^ = Scaling 
by ^ means 77 = 
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